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ABSTRACT 

Traditional pulsar polarization sweep analysis starts from the point dipole rotating vector model 
(RVM) approximation. If augmented by a measurement of the sweep phase shift, one obtains an 
estimate of the emission altitude (Blaskiewicz, Cordes, & Wasserman). However, a more realistic 
treatment of field line sweepback and finite altitude effects shows that this estimate breaks down 
at modest altitude ~ O.IRlc- Such radio emission altitudes turn out to be relevant to the young 
energetic and millisecond pulsars that dominate the 7-ray population. We quantify the breakdown 
height as a function of viewing geometry and provide simple fitting formulae that allow observers 
to correct RVM-based height estimates, preserving reasonable accuracy to R ~ 0.3Rlc- We discuss 
briefly other observables that can check and improve height estimates. 
Subject headings: methods: numerical — polarization — pulsars: general 



^ ■ 

o : 



> 

T-H 

CO 

T-H 

o 

CN 



1. INTRODUCTION 

After nearly a half century of pulsar observations, 
we still do not know the detailed location of the emis- 
sion zones in the neutron star magnetosphere. How- 
ever the general consensus is that the radio emission 
arises from the 'open' field line zone above the magnetic 
poles at modest altitudes, from a few to a few tens of 
neutron star radii. In contrast, the 7 -ray emission, as 
measured by Fermi (jAbdo et al.l l2010). is dominated by 
high altitudes > O.IRlc, where the light cylinder ra- 
dius is Rlc — cP/2tt. Thus the emission zones and 
light curves for these two bands generally differ. How- 
ever, recently Fermi has detected 7-ray emission from 
a number of millisecond pulsars where the entire mag- 
netosphere is outside of Rns ~ 0.2 /P ms Rl c so that 
radio emission m ust be from 'high altitude' ( Kerr et all 
1 20121). Further. iKarastergiou fc Johnston! (j2007[ ) and 
Johnston & Weisberg (|2QQ6f ) have found evidence that 
for young energetic pulsars, the radio emission is domi- 
nated by an altitude of - 1000 km (~ 100R NS ). This is 
~ 0.2Rlc for P=100 ms, and it is precisely such young, 
energetic pulsars which are 7-bright. Thus, if one is in- 
terested in 7-emitting pulsars, one must also consider 
radio emission from an appreciable fraction of the light 
cylinder radius. 

Since the first radio observations, the high linear po- 
larization and rapid position angle sweep of many pul- 
sars at cm wavelength have been used as a clue to the 
geometry of the emission zone. The fou n datio n for 
such study is the Radhakris hnan fc Cookej (|1969f ) 'ro- 
tating vector model' (RVM), which follows the sweep 
of the magnetic field line tangent of a point dipole as 
projected on the sky. Of course, finite altitude ra- 
dio emission violates the point source RVM a ssumption 
and lBlaskiewicz. Cordes fc Wasserman! (j 199 1\) (hereafter 
BCW) gave simple approximations for the effects of rel- 
ativist ic aberration at small altitude. In this approxima- 
tion, the polarization position angle is 
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where the inclination angle between rotation axis and 
magnetic axis is a, the viewing angle is £, and the pulse 
phase is (j). The RVM formula is recovered in the limit 
as the scaled emission height, r = r ern /RLC goes to 
zero. Here the principal effect is a lag in the phase of the 
maximum rate of the polarization sweep d^/d0| ma:c of 
Acj) « 2r from the phase of the magnetic axis. 

If the absolute position angle of the magnetic axis on 
the plane of the s ky is known (eg. from t he p osition angle 
of the spin axis), Hib schman fc Aronsl (j2001f ) show that 
the observed polarization gives a second height estimate, 
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(Dyks 2008). In practice it is generally unclear how to 
measure the magnetic axis polarization angles; most au- 
thors treat Aip as a nuisance parameter. 

Of course, both formulae presume knowledge of the 
phase of closest approach of the magnetic axis <j) = 0. 
The phase of the radio pulse peak is often used, but these 
pulses can have complex, multi-component morphology. 
Further, the special relativistic effects shift the intensity 
peak forward, giving a net observable lag of the polar- 
ization sweep from the intensity peak of Acj) w 4r. The 
shifts have been clarified and e xtended to include the ef - 
fect s of fie l d line sweepback bv lDvks fc Hardin g (2004), 
and iDvksl ([2008). Nevertheless, observers generally fit 
to the zero altitude (RVM) limits of the formula to con- 
strain a and £ and, when possible, estimate the shift 
of d^/d0| max to constrain the altitude, using the linear 
(BCW) scaling. While this works adequately for many 
non-recycled pulsars, relatively high altitude emission is 
inferred for young energetic objects. For millisecond pul- 
sars the basic RVM model often does not fit well. 

Thus, recent strong interest in 7-ray emitting pulsars 
draws our attention to objects where the radio emission 
may extend to O.IRlc or higher, where the standard 
RVM treatment is suspect. We seek here to quantify 
this breakdown: if one applies an RVM/BCW fit and 
obtains estimates of the magnetic inclination angle a/, 
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viewing angle and emission height r/, for what ranges 
of these parameters are these fits 'valid', i.e. when do the 
fit values and uncertainty ranges include (at some pre- 
scribed probability) the real value r r ? We develop this 
analysis as a guide to observers wishing to interpret pul- 
sar polarization data and as an indication to situations 
where detailed fits to numerical models (eg. Parent et al. 
2011) are required. In addition, we suggest analytic cor- 
rections to allow useful Tf estimates from simple RVM 
fits to extend to somewhat higher altitude. 

2. SIMULATION MODEL ASSUMPTIONS 

Our approach is to use a specific 3-D magnetosphere 
model with plausible radio emission zones, to 'fit' the re- 
sulting light curve and polarization sweep with the point 
dipole RVM formula and to parametrize the errors. For 
simplicity the field lines are given by the basic swept back 
(retarded ) dipole popular in models of high altitude 7-ray 
emission (Ro mani fc W atters 2010) and we assume that 
the radiating particle bunches follow the magnetic field 
lines. In the spirit of the RVM model, we make a simple 
geometric construction, projecting the field line tangent 
at the emission point in the lab frame onto the plane of 
the sky and assume that the radio emission is polarized 
parallel to (or perpendicular to) this vector. We do not 
attempt here to superpose multiple emission heights or 
to compute intrinsic polarization fractions. Nor do we 
include other physical effects such as possible cross-field 
drift of the emitting charge bunches, current-induced de- 
part ure from the vac uum structure for energetic pul- 
sars (Spitkovsky 2006)) or higher-order multipole/offset 
dipole effects that may be import ant in the small mag- 
netospheres of millisecond pulsars (Ha rding fc Muslimovl 
[20lTh . While our simple construction ignores these pos- 
sible effects, we do capture the dominant effect of dipole 
sweep-back and our computed polarization sweeps pass 
smoothly to the RVM model curves at low altitude; the 
other physical effects likely only dominate very close to 
the light cylinder. 

We assume here that the radio emission comes from 
a single altitude, within the open zone. We then must 
define the open zone shape and the illumination across 
it. Of course, there is a formal cap shape for the vacuum 
retarded dipole solution, where the locus of field lines 
tangent to the light cylinder trace to a cap on the sur- 
face with opening angle 0R(0 C ap) varying with azimuth 
(j) cap around the magnetic axis. Alternatively, it is com- 
mon to assume a simple circular cap, with surface angle 
^c(^cop) = constant. To roughly match the open zone 
beam sizes at an emission height of O.IRlc we chose 
a surface cap angle of 6c = 2° for a neutron star of 
Rns — 10 _3 i?LC 7 i-e. a ~ 0.2 s pulsar. 

For simplicity and to follow the BCW picture, we illu- 
minate the open zone with a simple Gaussian profile 

I oc e ( ^^ o)2 , with O = 2°/V / m~5 (3) 

so that the intensity falls by 5x at the 'edge' of the 
simple circular cap. The angles are measured at the 
star surface, although the corresponding radio flux may 
be emitted at high altitude. We note that there is 
some evidence that a conal intensity distribution with a 
patchy illumination may be more typical of many pulsars 
(iLvne and Manchester! 1 19881 : iKarastergiou fc Johnston] 
200% 



To generate a model polarization sweep we select a 
magnetic inclination, a r , and emission height, r r . We 
then project the swept-back field lines at this altitude 
to the plane of sky and record the results on a 2D sky 
map. Horizontal cuts across this map at a given viewing 
angle, £ r , give the polarization angle sweep, ip(4>). We 
assign 'measurement' errors to each value inversely pro- 
portional to the pulse flux at its phase. We assume that 
the observer's integration achieves a uniform signal-to- 
noise at pulse maximum, so that the polarization mea- 
surement error there is 1°. For pulsars observed far from 
the magnetic axis at large \/3\ = |£ — a\ this implies longer 
integration. As the pulse flux falls toward the edge of the 
open zone the polarization angle uncertainties increase. 

2.1. Estimating = 

Use of the simple Gaussian illumination with the pulse 
phase at the intensity peak (the projected phase of clos- 
est approach to the magnetic axis) corresponds to the 
BCW assumptions. Except for very high altitude emis- 
sion, where field lines overlap in the sky map and pulse 
caustics can occur, this gives a simple prescription from 
which = may be estimated via the BCW shift. How- 
ever, conal emission concentrated to the cap edge sig- 
nificantly complicates the determination of pulse phase. 
One effect is the variable sweep-back at the leading and 
trailing edge of the cap. Another is the particular shape 
of the open zone boundary. We illustrate these effects by 
marking a 'peak phase', the midpoint of the projected 
open zone boundary, both for a simple circular cap and 
for the more detailed retarded dipole cap. Figure [1] dis- 
plays the peak phase shifts for these different definitions. 
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Fig. 1. — Pulse phase estimates for a = 45°, </>(£) = for two 
different altitudes, r = O.IRlc and r = 0.3Rlc- Gc'- Peak phase 
from maximum of simple Gaussian intensity weighting. Gc'- Peak 
phase from center of cap edges (circular cap). Cr: Peak phase 
from center of cap edges (retarded dipole cap). 

Not unexpectedly, Figure [T] shows that the pulse phase 
is more sensitive to the details of the open zone geometry 
for a conal emission zone. The offsets shown there illus- 
trate the effect of the retarded potential field line flaring 
at high altitude. To this should be added the uncertain- 
ties associated with identifying the magnetic axis phase 
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in the presence of patchy conal emission and non-dipole 
field structure (for millisecond pulsars). Nevertheless, as 
we shall show, a substantial fraction of the phase off- 
set is insensitive to the choice of cap center, and can be 
corrected. 

2.2. 'Fitting' an RVM curve 

The retarded dipole field structure increasingly de- 
parts from the point dipole as the scaled emission height 
r — ^em/RLC approaches unity. Thus if one fits polar- 
ization data for a low altitude emitter with the RVM 
model, the fit parameters Q^/,C/? an d ff = at 
the x 2 minimum will be good approximations to the real 
values (a r ,C r ,r r ). For modest r r the RVM fit will ab- 
sorb the sweep shape departures, (correctably) biasing 
the parameter estimates, while retaining reasonable x 2 • 
At large altitude, the x 2 wm be poor, the parameters will 
be uncorrectably far from the true values and a fit to a 
detailed numerical model will be required. The key ques- 
tion is how, with realistic errors a pa, the unabsorbed 
distortion grows. 

We use our estimated a pa to construct a 'x 2 ' weighted 
departure of the RVM model from the detailed retarded 
field simulation. This is the weighted systematic error 
caused by the inability of the RVM model to absorb the 
detailed shape of the retarded field curve. In a real obser- 
vation, additional statistical measurement errors would 
increase 'x 2 ' above our model value, especially for small 
r. Any unmodeled physical effects should additionally 
increase the value of 'x 2 ' above ~ 1/ (degree of freedom) 
at the minimum. Observers typically adopt the increase 
Ax 2 = X 2 — Xmin t° estimate the confidence intervals 
on the fit parameters. We are free to do the same here, 
since our prescription weights appropriately show where 
the model parameters are most sensitive to the data val- 
ues. We have confirmed this by fits to a series of Monte 
Carlo simulations of polarization angle data with added 
statistical errors, showing that Ax 2 follows the usual dis- 
tribution for the appropriate numbers of degrees of free- 
dom. 
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CORRECTING FOR BIAS IN THE RVM HEIGHT 
ESTIMATES 



Our principal goal here is to test the utility of stan- 
dard RVM fits and to provide a prescription to allow 
these fits wider applicability for pulsars with high alti- 
tude emission. To do this we compare the RVM fit esti- 
mate 77 with the simulated value r r . Since the mapping 
is not simple, statements about ranges of validity are per- 
force statistical. This makes our answers mildly sensitive 
to the distribution in the underlying pulsar population. 
Here we assume that our parent pulsar population has 
isotropically distributed inclination and viewing angles, 
ie. Prob(a) ex sin(a), Prob(C) oc sin(£) while the al- 
titude is distributed uniformly on < r r < 0.3. Note 
that we only observe a usable polarization sweep if a 
pulsar produces a minimum number of phase bins (here 
A(j) bs > 0.1). In turn, this means that our observable 
pulsar population is biased toward modest \/3\ = \a — £|. 

We generate a set of pulsar models and apply 
the RVM fits. This delivers a set of observables 
C/j r /j °"q!/j °"C/' wnere the fit values are deter- 
mined by x 2 minimum and the error ranges are esti- 
mated from the curvature of the x 2 surface. An observer 



0.15 
o 0.1 
0.05 



— 

0.15 
o 0.1 
0.05 




— 




E B 




C c 



Go 




0.1 0.2 

Tf/RLC 



0.1 0.2 

Tf/RLC 



Fig. 2. — Altitude limits for effective RVM fits. Each panel 
shows the distribution of simulated model fits (color bar) in offset 
from the true altitude as a function of fit altitude Vf/R^c- The 
dark band shows the systematic bias in the fit offset. The four 
panels are for different assumptions about the cap illumination 
and method of estimating the true phase of <j> = 0. Eb- Perfect 
knowledge of the location of magnetic axis in phase, without the 
use of an intensity model. Gq'- Simple Gaussian intensity peak, 
<f> = inferred from the altitude dependent shift of Imax • Cc '• Peak 
intensity assigned to the center of a double pulse from edges of an 
open zone circular cap footpoints. Cr: Peak intensity assigned to 
the center of a double pulse from edges of an open zone above the 
detailed retarded dipole cap. The green curve shows our estimate 
of the bias, Equation (|4j). 



presented with this set of measurements must infer the 
original pulsar properties. 

Focusing here on the height measurement, we test the 
systematic bias in the RVM estimate. For best compari- 
son with the BCW assumptions, we work with the height 
determined from the phase lag measured from the peak 
of a Gaussian pulse centered on the magnetic axis. In 
Figure [2j the color scale represent the number of pulsars 
in the simulated population at a given altitude derived 
from fitting RVM versus the difference between fitted and 
real altitude. Figure [2] shows that 77 increasingly un- 
derestimates r r at increasing altitude. A simple formula 
to provide improved height estimates from RVM fits 
is then 

r' f =r f + 0.2(r f /0.5) 2 , (4) 

as plotted in Figure El The line fits best to the darkest 
ridge (The ridge that contains a majority of simulated 
pulsars) for the models using the maximum of a simple 
Gaussian intensity peak (Gc) and the center of the cap 
edges for a circular cap (Cc)- For the case using the cen- 
ter of the cap edge for a retarded dipole cap (Cr), the 
line slightly under-predicts the darkest ridge and does 
not capture the behavior of the second ridge which is 
caused by the shift of the central line from the cap notch 
(see Figure [1]). We can (unrealistically) assume that we 
know where in phase the magnetic axis is located and 
calculate the altitude from the shift in polarization di- 
rectly. Inaccuracies in altitude are then from Equation 
([2]) alone. With the assumption of perfect knowledge 
of the magnetic axis (Eb), we see the departure from 
the BCW formulation occurs at lower altitudes. Appar- 
ently, the estimate A</> = — 2 77 for the peak intensity 
shift preserves good accuracy to higher altitude than the 
Acj) = +2r/ shift of the polarization sweep, especially 



4 

when the intensity arises from a circular cap. 
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Fig. 3. — Maximum useful rf altitude (color bar) in the (af, C/) 
plane for four assumptions about the pulse intensity beam shape 
(see text for our criterion for good fit accuracy). Left: BCW esti- 
mates before correction. Right: corrected heights using Equation 
(UJ. Green contours indicate the area where at least fifty simulated 
pulsars were fit to an (ck/,C/) pair. 



In practice, the height offset depends on the geomet- 
rical angles (a, C)- I n addition, the height estimate is 
affected by uncertainty in estimating the polarization 
sweep lag, i.e. in determining the phase of the pulse 
(or equivalently the phase of the magnetic axis). These 
effects are shown in Figure 03 For each panel we show, 
as a function of the estimated angles (ctfXf), the maxi- 
mum height (color bar) at which the estimated altitude 
is accurate. For the estimate to be useful, we require 
that r r lies in the range rf ± cr rf for a large fraction 
(99%) of the observable model pulsars. At small altitude 
this is always true. At large altitude the distortion due 
to the retarded field structure causes increasing depar- 
ture from the BCW estimate. Once too small a fraction 
of models produce useful fits, the BCW approximation 
'breaks down'. Lowering the required fraction does not 
drastically change the results seen in Figure [3j since the 
fraction of failing models increases very rapidly with fit 
altitude. Also shown is a green contour that marks the 



area where the bins contain at least fifty simulated pul- 
sars. Uncolored bins are where the BCW approximation 
is inaccurate at the lowest altitude. The contours are 
independent of the intensity model (the contours are the 
same for each model) because the a/, Q bin depends 
only on the polarization sweep which is calculated inde- 
pendently of the intensity model. 

A strong dependence between the break-down altitude 
and af and Q exists as can be seen in Figure [U This 
is not due to any difficulty in finding the phase center 
but arises from the nontrivial relation between the shift 
in the maximum sweep of the polarization and the ge- 
ometry angles. In Figure [3j we can see that for af 
and (f further from 90°, BCW tends to break down at a 
lower altitude. The shift in the maximum sweep of the 
polarization angles for these values is smaller than pre- 
dicted by the BCW model. Since the BCW model has 
no dependence on a and £, it is not surprising that the 
break-down altitude has a dependence on these angles. 

The panels show the maximum useful height for four 
different estimates of the phase lag: (top-to-bottom) per- 
fect knowledge of the magnetic axis, a Gaussian pulse 
peaked on the magnetic axis field line, a 'conal' pulse 
from a field lines with a circular cap on the star and 
a 'conal' pulse with a cap determined by the detailed 
open zone of the retarded vacuum solution. Notice that 
most observed pulsars have modest \(3\ — \( — a|, and 
are close to the diagonal. The right panels show the 
equivalent maximum useful height when the estimate has 
been corrected according to Equation (j4j) . While the un- 
corrected estimates for the Gaussian pulse peak model 
are useful only to an average (over af and Q) height 
of fj = 0. 11 Rlc, the corrected estimates are usable to 
higher altitudes (reaching ~ 0.3 for the commonly 
observed case of near-orthogonal rotators) with an aver- 
age of YJ = 0.22Rlc- Again, corrected RVM estimates 
from a model radio pulse do better than estimates as- 
suming perfect knowledge of the magnetic axis, since the 
retarded potential phase shifts are a fractionally larger 
contribution to the phase offset in this case. 

We can improve the heuristic correction function by 
including the viewing geometry. The bulk of the sensi- 
tivity is evidently due to as illustrated by the rela- 
tively small dispersion of the rf error for individual Q 
slices (see Figure 2] for a Gaussian central pulse). Ac- 
cordingly, we have made an alternate corrected height 
estimate 

r f f = r f + [0.3 + 0.7| cos(()\](r f /0.5) 3 (5) 

where rf = A0/4, as usual. This greatly extends the 
range for which a simple RVM height estimate can be 
used (Figure [5j). This estimate, based on a Gaussian ra- 
dio pulse emitted along the swept back magnetic axis, is 
in general the best function for an observer to use with no 
other information. It provides significant improvement in 
the emission height accuracy for the circular cone pulse 
profiles. 

Of course if one has reason to believe that a particular 
pulse profile shape is more accurate, a different correction 
function may be preferred. For example if one had a 
double pulse arising from the open zone edges (Cr) and 
had high confidence that this pulse filled the retarded 
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Fig. 4. — Altitude limits for effective RVM fits. Each panel 
shows the distribution of simulated model fits (color bar) in off- 
set from the true altitude as a function of fit altitude r$ / Rlc for 
different assuming a simple Gaussian intensity peak, <j> = in- 
ferred from the altitude dependent shift of Imax- For these plots, 
the simulated pulsar population has been summed over ctf to em- 
phasize the dominance of £j in the correlation. The green curve 
shows our estimate of the bias with dependence on Equation 
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Fig. 5. — Maximum altitude for accurate height estimates (color 
bar) in the (a/, Cf ) plane, after applying Equation ((5} (see text for 
our criterion for good fit accuracy). Note that the improvement is 
best for a circular (Gaussian or conal) cap. Green contours indicate 
the area where at least fifty simulated pulsars were fit to an (ctf,Q) 
pair. 



(6) 



vacuum dipole open zone, one would correct by 

r' f =r f + [0.3 + 2|cos(C)| 2 ](r / /0.5) 3 . 

This formulation raises the average over af and Q of 
the maximum useful height from rj = 0.05Rlc with no 
correction to rj = 0.15Rlc- 



In general, we recommend that when an observer fits 
an RVM model to pulsar data, obtaining viewing angle 
and polarization sweep lag measurements, they correct 
their height estimate using Equation (j5j). This is partic- 
ularly useful whenever the RVM fit appears statistically 
adequate, but the resulting phase lag suggests a signifi- 
cant emission height. The change to the estimated height 
will be small for 77 < 0.2, but the accuracy of the result- 
ing estimate will be greatly increased. 

Of course, whenever x 2 /DoF ^> 1 at the fit minimum, 
it is a sign that the model is inadequate. In many cases, 
this will be due to unmodele d orthogonal mod e jumps 
and intervening scattering ([Karastergioul 120091 ), higher 
order multipoles, etc. However, for large altitudes and 
multi-altitude emission the effects of sweep back and the 
formation of caustics (which dominate 7-ray light curves) 
become dominant. The observer should be aware that 
large y 2 at the fit minimum can signal such effects and, 
when the inferred altitude is large, consider fitting the 
data to numerical models of 3-D pulsar magnetospheres. 

4. HEIGHT CALCULATION FROM SHIFT IN $ 

We can alternatively estimate r f a nd errors using the 
shift in ip (|Hibschman fc Aronsll2QQlf ). 
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r cos(a) 



+ - cos(C — a) 
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r sin(a) sin(C— a) 
(7) 



or, in the small \(3\ — \( — a\ limit, « cos(a). As 
before, we compute the residual, r r — 77, as a function of 
aj, (j, and Tf. To estimate an emission height from the 
polarization shift in 0, one needs an estimate for <fi — 0, 
e.g. from a pulse peak intensity model; no such intensity 
model is needed if we have a measurable shift in i/;. The 
increase with 77 are shown in Figure [6j where the left 
panel uses the small /3 limit while the right uses the full 
formula. As for the A</> estimate, the errors increase 
with rf. However here, even when the full equation ([7]) 
is used, the corrections show a substantial spread. In 
fact the uncorrected formula proves accurate (\rj — r r \ 
within <j rf 99% of the time) only for YJ < 0.08 (where YJ 
is again the average over af and Q) and for Q < 60° or 
Q > 120°. For near-orthogonal rotators the estimate is 
unreliable at the lowest altitudes. 




r f /R LC 



Fig. 6. — Altitude limits for effective RVM fits using the shift in 
Each panel shows the distribution of simulated model fits (color 
bar) in offset from the true altitude as a function of fit altitude 
rf /Rlc- The dark band shows the systematic bias in the fit offset. 
The residual is more scattered when the altitude is measured from 
the shift in i/j instead of the shift in cf) of the polarization sweep. On 
the left is the residual using the small \/3\ limit. The green curve 
shows our estimate of the bias from the shift in ip, Equation ([8]). 
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A heuristic correction to the Aip estimate for Equation 
can be made for Q < 60° or ( f > 120° 



r f +0A(r f /0.5)' 



(8) 



which allows accurate estimates to r/ = 0.12Rlc- In- 
cluding the £ dependence, 



r' f =r f + [0.2 + 0.1| cos(C)| 2 ](r / /0.5) 2 



(9) 



raises the useful range to rj = 0.18Rlc- Considering 
that the correction for the common orthogonal rotator 
case is especially poor, and that it is often difficult to 
infer the intrinsic ^o, height estimates from the phase 
shift remain much more useful. 

5. PULSE WIDTH DEPENDENCE ON EMISSION HEIGHT 

Since the field lines flare in the open zone, the full phase 
width W of the observed radio pulse can also be checked 
against the expected radio emission altitude. The stan- 
dard prescription assumes a circular cap and static dipole 
field lines to infer a minimum height 
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- arccos 



cos(a) cos(C) 



sm(a) sm(C) cos [ — 



(10) 

In Figure [7] we show that the retarded dipole field flares 
more than predicted by this simple formula and hence the 
minimum height in Equation ([TO]) is an over-estimate. 
Thus, in general, lower altitudes are consistent with a 
given observed pulse width than suggested by this for- 
mula. Moreover, we expect that the general effect of 
currents in the magnetosphere will be to increase the 
foot-point angles of the open zone. This further increases 
the allowed If at a given height. 
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Fig. 7. — Altitude limits for effective pulse width. Each panel 
shows the distribution of simulated model fits (color bar) in offset 
from the true altitude as a function of fit altitude rf /Rlc- Wc'- 
Circular cap. Wr: Retarded dipole cap. The simple static dipole 
formula overestimates the altitude needed to accommodate a given 
pulse with W in the open zone. The error depends on the viewing 
geometry a and £, and bifurcates for the 'notched' cap of the formal 
retarded potential open zone. 



In general, larger widths are still most easily accom- 
modated at large r or small a, but sweep-back and mag- 
netospheric currents substantially weaken the minimum 
altitude constraints from the commonly used Equation 
([TO]) . Given the large sensitivity to the details of the 
open zone volume and the presently unknown effect of 
magnetospheric currents, it is not worth developing cor- 
rections to this formula. 



6. CONCLUSIONS; EXAMPLES FROM LITERATURE 

We conclude by examining a few RVM/BCW estimates 
of e mission height present in the literature. 

In iRomani et al.l ([20111 ), estimates were used to 
suggest large emission heights for two young energetic 
pulsars. For J0538+2817 the shift gives r f = 0.15R LC , 
but RVM fitting only weakly constraints £. Applying 
Equation (J2J), we would infer = 0.17 Rlc, a small, 
but significant increase which makes it easier to accom- 
modate the large observed pulse width. Similarly PSR 
J1740+1000 gives rf = 0.12Rlc- Here we constrain 
C = 80° to 130°, so that the corrected fit altitude (Equa- 
tion [4] or Equation [5]) is = 0.13.RLC, again a small but 
statistically significant increase. 

For mi llisecond pulsars t he effects can be larger. For 
example. [Keith et al.l([2012[ ) find that RVM fitting can be 
usefully applied to several recycled pulsars. PSR J1502- 
6752 (P=26.7ms) is a mildly recycled pulsar for which 
the phase lag implies rf = 0.2 Rlc- With no signifi- 
cant ( constraints, we apply Equation ([4]) to infer a 16% 
altitude increase to r' f = 0.23R L c- Similarly PSR J1708- 
3506 (P = 4.5 ms) has a phase shift implying 0.19Rlc, 
which we correct to 0.21 — 0.22Rlc- For this pulsar, a 
naive application of the pulse width formula ([TO]) gives 
altitudes of rw 10 — 0.65Rlc (10% peak width). How- 
ever, the increased and decreased pulse width height 
from sweepback effects (Figure [7j), along with additional 
current-induced open zone growth, make it likely that 
the pulse width can be accommodated at the corrected 
h eight. 

iKeith et al.l ([20121 ) also report a RVM/BCW height 
r f = 0.UR LC for the P=2.7ms pulsar J1811-2404, along 
with well constrained viewing angles of a = 89.7° and 
j3 = 21°. While our full analysis does not cover this alti- 
tude, as Figure [5] shows the corrections of Equation (|5j) 
give a very high accuracy for orthogonal rotators viewed 
near 90°. Note in Figured! bottom panel, that the cor- 
rection function is nearly linear thus extrapolation to 
somewhat higher values may be justified. Naively apply- 
ing this correction we get = 0.81 Rlc- We certainly 
cannot trust this value in detail since plasma effects and 
other perturbations may be relevant at such altitudes. 
However, the correction is certainly large and it brings 
the expected height up to an altitude where the very 
wide observed radio pulse, and the likely detection of 
emission from both open zones, can be easily accommo- 
dated. Certainly simple RVM/BCW fitting is inadequate 
for this pulsar and one should use a detailed model for 
the high altitude field geometry. 



Our exercise extends the range of utility of RVM-fit 
polarization sweeps for inferring the altitudes of radio 
pulsar emission. For fit altitudes less than rf — 0.25Rlc 
the corrections are not large, but they are systematic 
and, for high S/N data localizing the phase of maximum 
polarization sweep, they can be highly significant. We 
thus believe it is worth applying our recommended cor- 
rection. For larger altitudes the corrections grow rapidly, 
but we caution that as one approaches the light cylinder, 
current-induced distortions should increase and, except 
for near-orthogonal rotators, one would expect the RVM 
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formulae to provide a poor fit in any case. Fitting to 
detailed numerical models is then preferred. In all cases 
the dominant residual uncertainty is likely in locating 
the phase of the radio pulse. We also checked the use 
of absolute polarization axis position angles and pulse 
width to constrain the emission height. Here the difficul- 
ties in establishing the unperturbed ipo and the expected 
distortions of the open zone boundaries by currents, etc. 
make the estimates much less useful. Nevertheless, we 
have shown that the effects of sweepback do go in the 
direction of reconciling observed pulsar properties to a 
consistent emission height: larger heights are inferred by 
a given Aif; shift and larger pulse widths can be accom- 



modated at a given height. We feel, however, that the 
corrections are less quantitative than for A<fr. 

In sum, since observers will continue to apply analytic 
RVM fits to pulsar polarization data, by applying our 
recommended correction (Equation [5j), these results can 
continue to give accurate height estimates to < 0.3Rlc- 
At higher heights which will be common for millisecond 
pulsars, a fit to more detailed numerical models is likely 
warranted. 



This work was supported in part by NASA grants 
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